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We give a summary about the various contributions which have to be calculated 
in order to obtain the next-to-leading logarithmic result for the branching ratio 
BR{B Xs-y). Combining all these ingredients, which were obtained by dif- 
ferent groups, a complete next-to-leading-logarithmic prediction of the inclusive 
decay rate was recently presented in the literature. The theoretical uncertainty 
in the partonic decay rate is now at the 10% level, i.e., less than half of the error 
in the previous leading-logarithmic result. We also mention the impact of non- 
perturbative corrections which scale like 1 / and discuss in some more detail the 
recently discovered corrections which scale like It turns out that the 

and the terms lead to corrections to the branching ratio BR{B — ► Xs^) well 

below the 10% level. 



1 Introduction 

The B Xs'y decay has found increasing attention over the last ten years. It 
provides an alternative approach in the search for physics beyond the standard 
model (SM). This decay, like other rare B meson decays, does not arise at 
the tree-level in the SM but is induced by one-loop W-exchange diagrams, so 
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nonstandard contributions (charged scalar exchanges, SUSY one-loop diagrams 
etc.) are not suppressed by an extra factor a/An relative to the standard model 
amplitude. This high sensitivity for nonstandard contributions implies the 
possibility for an indirect observation of new physics, a strategy complementary 
to the direct production of new particles. The B X^^f decay plays already a 
very important role in restricting the parameter space of exlensions of the SM 
like the minimal supersymmetric standard model (MSSM)El'El. However, even 
within the SM, the B Xg'y decay is important for constraining the Cabibbo- 
Kobayashi-Maskawa matrix elements involving the top-quark, in particular 
\Vts\. For both reasons, precise experimental and theoretical work on this 
decay mode is required. 

On the theoretical side, the accuracy in the dominating peiitndiatiye con- 
tribution was recently improved to next-to-leading precision BqIj'BlHj: The 
renormalization scale dependence of the previous leading-log result at the 
±25%-level was substantially reduced to ±6% and the central value was shifted 
out-side the la bound of the CLEO measurement. Furthermore, the analysis 
of nonperturbtive contributions to the B — > Xs"f decay mode was also recently 
improved: The inclusive B —> Xs^ mode is theoretically much cleaner than 
the corresponding exclusive channels because no specific model is needed to 
describe the final hadronic state. According to Heavy Quark Effective The- 
ory the class of non-perturbative effects which scales like is expected 
to be well below 10% □. This numerical stalpmpaJJinlds also for the recently 
discovered non-perturbative contributions E30't3'E3E3 which scale like 
Thus the inclusive B Xg'y mode is well approximated by the partonic decay 
rate T{h — > Xg^) which can be analyzed in renormalization group improved 
perturbation theory. 

Before reporting on these theoretical improvements in detail, we summa- 
rize the experimental status: The observation of the exclusive B — )■ K*^ mode 
by CLEOEj in 1993 was the first evidence for a penguin decay ever. An updated 
value0 for the branching ratio is BR{B K*j) = (4.2 ± 0.8 ± 0.6) x 10"^ 
In 1994 the CLEO collaboration measured the inclusive B Xsj branching 
ratio to be (2.32 ± 0.57 ±^0.35) x lO"* where the first error is statistical and 
the second is systematic Ej: There are two separate CLEO analyses. The first 
one measures the inclusive photon spectrum from B-decay near the end point. 
The second technique constructs the inclusive rate by summing up the possi- 
ble exclusive final states. The branching ratio stated above is the average of 
the two measurements, taking into account the correlation between the two 
techniques. 

Thpe is also data from the LEP experiments: While DELPHI El in 1996 
and L3lla in 1993 have published the upper bounds BR{b sj) < 5.4 x 10"** 
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and BR(h —>■ 57) < 1.2 x 10 respectively, the preliminary measurement 
BR{b 57) = (3.38 ± 0.p.± 0.85) x lO"'^ by the ALEPH group was reported 
in the talk by F. ParodiEj at the_1997 Moriond meeting. A similar number 
was also quoted by T. SkwarnickiEil in the heavy flavor meeting held in Santa 
Barbara in July 1997. 

More precise measurements are expected from the upgraded CLEO detec- 
tor, as well as from the B-factories presently under construction at SLAC and 
KEK. In view of the expected high luminosity of the B-factories, experimental 
accuracy of below 10% appears to be in reach. 

The rest of the paper is organized as follows: Section 2 is devoted to the 
partonic (=perturbative) contribution to BR{B — > Xgj). We explain in some 
detail the various calculational steps leading to the next-to-leading logarithmic 
result. In section 3 we briefly discuss the impact of the recently discovered non- 
perturbative corrections which scale like 1/m^. 

2 Next-to-leading logarithmic corrections for B Xgj 

It is well-known that the QCD corrections enhance the partonic decay rate 
r(6 — > 57) by more than a factor of two. These QCD effects can be attributed 
to logarithms of the form a"(mf,) log™ (mb/M), where M = m* or M = mw 
and m < n (with n = 0, 1, 2, ...). In order to get a reasonable result at all, one 
has to sum at least the leading-log (LL) series {m — n). Working to next- to- 
leading-log (NLL) precision means that one is also rcsumming all the terms of 
the form as{mb) {a^imb) ln"(mb/M)). 

An appropriate framework to achieve the necessary resummations is an 
effective low-energy theory, obtained by integrating out the heavy particles 
which in the SM are the top quark and the W^-boson. The effective Hamiltonian 
relevant for 6 — s- 37 and 6 — s- in the SM and many of its extensions reads 

i?e//(6-s7) = -^A,^a(/i)0.(M) , (1) 

where Oi{fj,) are the relevant operators, Ci{fi) are the corresponding Wilson 
coefficients, which contain the complete top- and W- mass dependence, and 
At = VtfcVj* with Vij being the CKM matrix elements Neglecting operators 
with dimension > 6 which are suppressed by higher powers of l/m^/^ and 
using the equations of motion pEor the operators, one arrives at the following 
basis of dimension 6 operators c2l 

'^The CKM dependence globally factorizes, because we work in the approximation \u = 0. 
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0-2 = {cLaYhLa) {sL/Sj^iCLp) , 

O7 = {e/l67r'^)sacr'^'''^ imb{^l)R + ms{^i)L)ba Ff^^ , 

Os = {gs/16TT^)s^<j''''{mt{ti)R + ms{ti)L){X^^/2)bf,G^, . (2) 

Because the Wilson coefRcients of the penguin induced four-fermion operators 
O3, Og are very small, we do not list them here. In this framework the next- 
to-leading logarithmic terms as{mb) (a"(m6) log"'(mf,/mvK/t)) the b ^ 
amplitude have two sources: 

1 • The NLL Wilson coefRcients Ci{fi) at the scale /i w mi, contain leading 
and next-to-leading logarithmic terms in resummed form. 

2 • The 0{as) corrections to the matrix elements of the operators Oi yield 
next-to-leading order terms when multiplied by the (leading logarithmic part 
of the) Wilson coefRcients. 

We stress that only the sum of these two sources is independent of the renormal- 
ization scheme. Let us discuss in some more detail the contributions mentioned 
in 1 and 2: 

ad 1 • From the /i-independence of the effective Hamiltonian, one can 
derive a renormalization group equation (RGE) for the Wilson coefRcients 

tJ'-^Ciin) = jj^Cjifi) , (3) 

where the (8 x 8) matrix 7 is the anomalous dimension matrix of the operators 
Oi. To solve this first order differential equation one explicitly needs initial 
conditions Cilpo) at some scale /io as well as the anomalous dimension matrix 

la: The initial conditions are obtained by matching the effective theory to the 
full standard model theory at the scale fiQ = where denotes a scale 
of order mw or mt- At this scale, the matrix elements of the operators in 
the effective theory lead to the same logarithms as the full theory calculation. 
Consequently, the Wilson coefficients Ci{^w) only pick up small QCD correc- 
tions, which can be calculated in fixed-order perturbation theory. In the LL 
(NLL) program, the matching has to be worked out to order (a^) precision, 
lb: Solving the RGE (|) and using the Ci{^w) of Step la as initial condi- 
tions, one performs the evolution of these Wilson coefRcients from /i = 
down to fj, = fib, where fib is of the order of mb. As the matrix elements of 
the operators evaluated at the low scale fib are free of large logarithms, the 
latter are contained in resummed form in the Wilson coefficients. For a LL 
(NLL) calculation, this RGE step has to be performed using the anomalous 
dimension matrix 7^ up to order a] {aj.). 
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(a) (6) (c) 



Figure 1: a) Typical diagram (full theory) contributing in the NLL matching calculation, 
b) Typical diagram contributing to the matrix element of the operator 02- c) Typical 
contribution to the 0{a^) anomalous dimension matrix. 

ad 2 • The matrix elements of the operators {sj\Oi{^)\b) at the scale ^ = fib 
have to be calculated to order (al) in the LL (NLL) calculation. 

Until recently, onlx the leading logarithmic (LL) perturbative QCD were 
known systematically H The error in this approximation was dominated by a 
large renormalization scale-dependence at the ±25% level. The measurement 
of the CLEO coUaborationU overlaps with the estimates based on leading log- 
arithmic calculations (or with some next-to-leading effects par tially included) 
and the experimental and theoretical errors are comparableB'EjOuj. However, 
in view of the expected increase in the experimental precision in the near fu- 
ture, it became clear that a systematic inclusion of the NLL corrections was 
necessary. This ambitious NLL enterprise was recently completed. All three 
steps (la,lb,2) involve rather difficult calculations. The most difficult part 
in Step la is the two- loop (or order as) matching of the dipole operators O7 
and Os- It involves two- loop diagrams both in the full and in the effective 
JJieory (see Fig. |^a). This matching calculation was done by Adel and Yao 
□ some-time ago. As this is a crucial step in the NLL program, Greub and 
Hurth 13 recently confirmed their findings in a detailed re-calculation, using 

a somewhat different method. In order to match dimension 6 operators O7 

2 

and Os, it is sufficient to extract the terms of order mf, ^ (M = tovi/,^'*) 
from the standard model matrix elements for 6 — > 57 and b — > sg. Terms 
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supressed by additional powers of mb /M correspond to higher dimensional op- 
erators in the effective theory. In El the finite parts of the two-loop diagrams 
in the SM were calculated by means of the well-known Heavy Mass Expansion 
(HME) which naturally leads to a systematic expansion of Feynman diagrams 
in inverse powers of M. We mention here that the evolution of the Wilson 
coefficients between fi ~ rrit and ji = m^, to LL precision implied an addi- 
nal contribution of +10% to the leading- log prediction for the decay rate 
Most of this iicmtribution is automatically included in tli£_|NLL matching 
at the mvy-scaleCTD, because the first term of the LL-sum ofE3 is reproduced 
and higher order terms (a^ log(^^i^))"(n > 1) are rather small. In addition, 
the NLL matching result includes the first term of the NLL-sum. 

Step 2 basically consists of Bremsstrahlung corrections and virtual cor- 
rections. While the Bremsstrahlung corrections (together with some virtual 
corrections needed to cancel infrared singularities) were worked out some tipie 
ago by Ali and Greub a and have been confirmed and extended by Pott B, a 
complete analysis of the virtual corrections (up to the contributions of the 
four-fermion operators with very small coefhcients) was presented by Greub, 
Hurth and Wylercl. This calculation also involves two- loop diagrams where 
the full charm quark mass dependence has to be taken into account. A typical 
diagram is shown in Fig. By using Mellin-Barnes techniques in the Feyn- 
man parameter integrals, the result of these two-loop diagrams was obtained 
in the form 

/ 2 \ " 2 

C0+ E c„„ Pi log"^ , (4) 

n=0,l,2,...;m=0,l,2,3 ^"^b J "^b 

where the quantities cq and c„m are independent of mc- Note, that a finite 
result is obtained in the limit uric 0, as there is no naked logarithm of 
rn^/rn^. This observation is of some importance in the b ^ dj process, where 
the u-quark propagation in the loop is not CKM suppressed. It is, however, 
even more important that the inclusion of the 0{as) matrix elements leads to 
a drastic reduction of the renormalization scale uncertainty from about ±25% 
to about ±6%. Analytically, the reason is, that the term aslog(/i/TOb) which 
dominates the /i-dependence of the LL result, is cancelled by a corresponding 
term appearing in the 0{as) matrix element. Finally, the anomalous dimension 
matrix, (at 0{a'^)), Step lb, has been worked out by Chetyrkin, Misiak and 
MiinzQ. The calculation of the elements 7i7 and 7^8 (« = 1, ■■•,6) in the 0{a^) 
anomalous dimension matrix involves a huge number of three loop-diagrams 
from which the pole parts (in the d — 4 expansion) have to be extracted. For 
a typical diagram see Fig. Wc. The extraction of the pole parts were simplified 
by a clever decomposition of the scalar propagator. Moreover, the number 
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of necessary evanescent operators were reduced by a new choice of a basis of 
dimension 6 operators. Using the matching result (Step la), these authors 
obtained the next-to- leading correction to the Wilson coefficient CT{^b) which 
is the only relevant one for the b — > Xs"i decay rate. Numerically, the LL 
and the NLL values for CT{^b) are rather similar; the NLL corrections to the 
Wilson coefficient CrCub) lead to a change of the b Xg^ decay rate which 
does not exceed ±6%Ll: The new contributions can be split into a part which 
is due to the order corrections to the matching (Step la) and into a part 
stemming from the improved anomalous dimension matrix (Step lb). While 
individually these two parts are not so small (in the NDR scheme, which was 
used inEI), they almost cancel when combined as illustrated inQ. This shows 
that all the three different pieces, la, lb, 2, are numerically equally important. 

Combining the NLL calculations of all the three steps (la+b,2), the first 
complete theoretical prediction to NLL pecision for the 6 — > + 7 branching 
ratio was presented ini BR{B Xsj) = (3.28±0.33) x 10"'*. The error is due 
to the ±6% renormalization scale uncertainty and due to the ±8% combined 
uncertainty in the input parameters. 

3 1/m^ and 1/m^ corrections 

Neglecting perturbative QCD corrections and assuming that B Xs^ is due 
to the operator O7 only, the calculation of the differential decay rate basically 
amounts to work out the imaginary part of the forward scattering amplitude 

nq) 

T{q) J d'^x{B\TO:^{x)O7{0)\B) exp{iqx) . (5) 

Using the operator product expansion for TOj {x) Oj{0) and Heavy Quark 
Effective Theory methods, the decay width T{B ^s^) reads □ (modulo 
higher terms in the 1/mt, expansion) 

(07,07) _ IT/ 12^2. N A , 

SZ^, = 1^1-1^2 , (6) 

where Ai and A2 are the kinetic energy- and the chromomagnetic energy pa- 
rameters. Using Ai = ~0.5GeV^ and A2 = 0.12 GeF^ one gets (5^^ ~ -4%. 
As also the semileptonic decay width gets 1 / corrections which are negative 
(see e.g. E2I), these non-perturbative corrections tend to cancel in the branch- 
ing ratio BR{B — s- Xs^) and only about 1% remains. This contribution was 
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(o) (6) 



Figure 2: a)Feynman diagram from which the operator O arises, b) Relevant cut-diagram 
for the (O2, 07)-interference. 



already included in the theoretical NLL prediction presented in section 2 of 
this article. 

Recently, Voloshin0 considered the non-perturbative effects when includ- 
ing also the operator 02- This effect is generated from the diagram in Fig. 
^ (and from the one not shown where the gluon and the photon are inter- 
changed); 5 is a soft gluon interacting with the charm quarks in the loop. Up 
to a characteristic Lorentz structure, this loop is given by the integral 

\^ f'-\y ^ ^ 

Jo ~ kgx[l ~ x) — 2xykgkj 

As the gluon is soft, i.e., kg,kgkj ~ A^'^'^ mh/2 <C m^, the integral can be 
expanded in kg. The (formally) leading operator, denoted by O, is 

O - ^V,bV;,C, 57^(1 - j,)gsG,xb e>^^'^d^F,, . (8) 

Working out then the cut diagram shown in Fig. ^jb, one obtains the non- 
perturbative contribution F^^^-*^ to the decay width, which is due to the 
(O2, O7) interference. Normalizing this contribution by the LL partonic width, 
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one obtains 



_2_^^_^^^,,+0.03 . (9) 

Including this correction with the sign found inEl, the NLL prediction for the 
branching ratio becomes BR{B X^j) = (3.38 ± 0.33) x 10 
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As the expansion parameter is mbP^Qc d / fric ~ 0-6 (rather than ^qcdI 
it is not a priori clear whether formally higher order terms in the r n.^ expan- 
sion are numerically suppressed. More detailed investigationsE3E3o show that 
higher order terms are indeed suppressed, because the corresponding expansion 
coefficients are small. 

We mention that the analogous 1 / ral effect has been found independently 
in the exclusive mode B K*^ in ref. EJ. Numerically, the effect there is also 
at the few percent level. 



4 Summary 

Collecting all NLL contributions and the small nonperturbative correction 
which scales with m^, the final analysis done by Chetyrkin, Misiak and Miinz 
yields BR{B X^-/) = (3.38 ± 0.33) x 10"'' when also the +3% shift due to 
the non-perturbative effects from the l/m^ corrections is included. The theo- 
retical error in the NLL prediction is reduced by a factor of 2 when compared 
with the LL result. This theoretical value for the branching ratio is in agree- 
ment with the CLEO measurement (at the 2cr-level) and also with the recent 
(preliminary) measurement by ALEPH. Clearly, the inclusive B — > Xg + 7 
mode will provide an interesting test of the SM and its extensions as soon as 
more precise experimental data become available. 



Note added: When fini^ing this article, we received the new work by 
Buras, KwiatkowsM and PottCJ. While these authors fully confirm the match- 
ing conditions byQa, their analysis is slighly different, leading to the branching 
ratio BR{B -> ^^7) = (3.48 ± 0.31) x 10"'^. The shift in the central value 
is due to systeroatically discarding next-next-leading order terms, while in the 
earlier analysisQLl some terms of this order were included. Also their estimate 
for the remaining renomalization scale dependence is somewhat different: The 
/i-uncertainties in the decay width for the radiative decay and the semileptonic 
decay were treated independently and added in quadrature. In the old analysis 
eH3 the scale /j. was varied simultaneously in both decays. As the semileptonic 
decay width is an increasing function of ^ while the radiative decay width is 
decreasing, a larger /i-uncertainty was obtained which as a more conservative 
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estimate we finally prefer. The results are fully compatible after all. 
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